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Abstract 

A new set of supersymmetric benchmark scenarios has recently been proposed in the 
context of the constrained MSSM (CMSSM) with universal soft supersymmetry-breaking 
masses, taking into account the constraints from LEP, b s'-f and — 2. These points have 
previously been used to discuss the physics reaches of different accelerators. In this paper, 
we discuss the prospects for discovering supersymmetric dark matter in these scenarios. We 
consider direct detection through spin-independent and spin-dependent nuclear scattering, 
as well as indirect detection through relic annihilations to neutrinos, photons, and positrons. 
We find that several of the benchmark scenarios offer good prospects for direct detection via 
spin-independent nuclear scattering and indirect detection via muons produced by neutrinos 
from relic annihilations inside the Sun, and some models offer good prospects for detecting 
photons from relic annihilations in the galactic centre. 
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1 Introduction 



After the closure of LEP, at the start of Run II of the Tevatron Collider, with the LHC 
experimental programme being prepared, and linear e^e~ collider projects under active 
discussion, now is an appropriate time to review the available experimental constraints on 
supersymmetry and assess the prospects for its discovery. In parallel with present and future 
accelerator projects, many non-accelerator experiments that may contribute to the search 
for supersymmetry are underway or in preparation. These include direct searches for the 
elastic scattering of astrophysical cold dark matter particles on target nuclei, and indirect 
searches for particles produced by the annihilations of supersymmetric relic particles inside 
the Sun or Earth, in the galactic centre or in the galactic halo. 

A set of benchmark supersymmetric model parameter choices was recently proposed IQ 
with the idea of exploring the possible phenomenological signatures in different classes of 
experiments in a systematic way. The proposed benchmark points were chosen by first 
implementing the constraints on the CMSSM parameter space imposed by previous ex- 
periments, such as the searches for sparticles and Higgs bosons at LEP and elsewhere, 
the measured rate for 6 — > 57 decay , and (optionally) the value oi — 2 recently reported 
by the BNL E821 experiment [^]. The CMSSM parameter space was also constrained by re- 
quiring the calculated supersymmetric relic density to fall within the range 0.1 < Vt^h"^ < 0.3 
preferred by astrophysics and cosmology. Four general regions of allowed parameter space 
were identified: a 'bulk' region at relatively low itlq and mi/2, a 'focus point' region |0, H at 



relatively large mo, a coannihilation 'tail' extending out to relatively large mi/2 0, [lOl' ^^'^ ^ 



possible 'funnel' between the focus point and coannihilation regions due to rapid annihilation 



via direct-channel Higgs boson poles W^ - 

The benchmark points were chosen not to provide an unbiased statistical sample of the 
CMSSM parameter space, which is in any case difficult to define in the absence of any unbi- 
ased a priori measure, but rather to select representative examples of different possibilities 
that cannot yet be logically excluded. Note that while these scenarios are confined to the 
context of supergravity, they span a large range of dark matter properties. While other 
supersymmetry-breaking schemes lead to a variety of collider signals, with respect to dark 
matter, they often predict vanishing or highly suppressed thermal relic densities for the most 
natural candidate, the neutralino. These alternative scenarios therefore typically have no 
viable dark matter candidates, at least without additional structure and an accompanying 
loss of predictability. 

Of the 13 benchmark points, B, C, G, I, and L lie within the 'bulk' region; E and F are 
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Figure 1: Benchmark points in the (a) {mQ,mi/2) and (b) (|/i|,Mi) planes. 



in the focus point region; A, D, H, and J are strung out along the coannihilation tail; and K 
and M are chosen at (relatively) large mi/2 and mo, in the rapid annihilation funnel regions. 
About half of the proposed points yield a value oi — 2 within two standard deviations of 
the value reported by BNL E821, but we did not impose this as an absolute requirement. For 
example, two points with < 0, the sign disfavoured by — were retained. Fig. ^provides 
an overview of the locations of the benchmark points in the (mo, mi/2) and (|/i|. Mi) planes. 
We see that the proposed scenarios mainly have mi/2 > fno, except for the two focus point 
models E and F. These also have larger values of Mi/|/i|, and therefore more Higgsino-like 
lightest supersymmetric particles (LSPs). Table displays many properties of the proposed 
scenarios, including the LSP mass, its gaugino composition, its cosmological relic density, 
and rates for the many astrophysical signatures to be discussed in subsequent sections of 
this paper. 

It was found previously [|l| that, in the (7^ — 2-friendly scenarios, supersymmetry was 
relatively easy to discover and study at future colliders such as the LHC and a linear collider 
with EcM = 1 TeV, which would be able to observe rather complementary subsets of CMSSM 
particles. However, some of the other points might escape detection, except via observations 
of the lightest neutral Higgs boson of the CMSSM. The most difficult points were typically 
those in the focus point region, at the tip of the coannihilation tail, or along the rapid- 
annihilation funnels, with points F, H, and M being particularly elusive. 

In this paper, we report on the prospects for the direct and indirect detection of as- 
trophysical dark matter for each of these benchmark points. We present cross sections for 
the elastic scattering of supersymmetric relic particles off both protons and neutrons via 
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Table 1: Parameters and dark matter observables for the benchmark points. The super- 



symmetric mass spectra are obtained using ISASUGRA 7.51 with the listed input param- 
eters. For all the benchmark points, we assume Aq = and rrit = 175 GeV. All masses 
and energies are in GeV. We define the gaugino fraction of the lightest neutralino x 



Rx = 



■'XI I 



+ |Z 



X2l 



^, where x = ^xi^ + ^X2^° + ^xs^u + ^xi^d- The neutralino relic 
density fi^^'^ taken from Table 2 of and were calculated using SSARD j7^/ . The spin- 
independent (spin- dependent) cross sections on protons ap (cr^p) and neutrons ((^n) ^'^^ 



muon 



calculated with Neutdriver / [7^ and are given in units of 10 pb (10 ph). The 
fluxes from the Sun (^®) and the Earth (^'^) are in units o/ km~^ yr""*^- The integrated 
photon fluxes (^^) for photon energy threshold Eth = 1 GeV (Eth = 50 GeVj are in 
units of 10~^^ cm~^ (10~^^ cm~^ ^~^)- Finally, S/B is the maximal value of the positron 
signal-to- background ratio, and Eopt is the energy at which this value is realized. 
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both spin-independent and spin-dependent matrix elements, the rates for observing muons 
induced by the colhsions in rock of energetic neutrinos produced by rehc annihilations inside 
the Sun and Earth, the rates for photons produced by annihilations in the galactic centre, 
and the rates for positrons produced by the annihilations of relic particles in the galactic 
halo. In all cases, we take into account the sensitivities of present and planned detectors 
in estimating the observability of signals from relic particles. We emphasize that all our 
results necessarily depend on the halo model used: this is particularly true for the photon 
signal from the galactic centre. This model-dependence enters when comparing the power of 
various experimental probes. However, for any given signature, our conclusions concerning 
the relative ease with which different models can be seen should be quite reliable. 

The structure of this paper is as follows. In Sec. ^ we review briefly the experimen- 
tal constraints that were used as inputs when proposing the benchmark points studied in 



this paper. In Sec. ^ we compare the predictions of two different codes, Neutdriver |14| 
and SSARD for direct dark matter detection, obtaining very similar results. We use 
Neutdriver to calculate muon rates from the Sun and Earth in Sec. ^, and we follow the 
analysis of to determine the photon and positron rates in Sees. ^ and ^, respectively. 
Finally, in Sec. |^ we draw some tentative conclusions about the detectability of dark matter 
particles in the different allowed regions of parameter space, and we contrast the prospects 
in accelerator and non-accelerator experiments. 



2 Constraints used to Select Benchmark Points 

We restrict our attention to a constrained version of the MSSM (CMSSM) which incor- 
porates a minimal supergravity-inspired model of soft supersymmetry breaking. Universal 
gaugino masses mi/2, scalar masses mo (including those of the Higgs multiplets) and tri- 
linear supersymmetry breaking parameters Ao are used as inputs at the supersymmetric 
grand unification scale. In this framework, the Higgs mixing parameter /i can be derived 
(up to a sign) from the other MSSM parameters by imposing the electroweak vacuum con- 
ditions for any given value of tan f3. Thus, given the set of input parameters determined by 
{mi/2, mo, Aq, tan (3, sgn(yu)}, the entire spectrum of sparticles can be derived. For simplicity, 
we further restrict our attention to Aq = 0. 

The available experimental and phenomenological constraints on the CMSSM parameter 
space were implemented in These include the experimental constraints obtained from 
searches for sparticles ||^ and Higgs bosons at LEP In particular, attention was restricted 
to parameter choices which guaranteed chargino masses m^± > 103.5 GeV and selectron 
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masses mg > 99.4 GeV |T^. The lower limit on the mass of a Standard Model Higgs boson 
imposed by the combined LEP experiments is 113.5 GeV and this limit also applies to 
the lightest supersymmetric Higgs boson h in the CMSSM. To calculate nih theoretically, we 
use the FeynHiggs code |TB[, which includes one-loop effects and also the leading two-loop 
contributions. To account for uncertainties in theoretical calculations of rrih [|TB|, for any given 
value of rrit, we restrict our CMSSM parameter choices to those yielding nih > 113 GeV. 
In addition, the theoretical value of nih in the MSSM is quite sensitive to m^, the pole 
mass of the top quark: we use = 175 GeV as default. All but one of the benchmark 
points satisfy nih > 113 GeV. In view of the expected accuracy ~ 3 GeV of the FeynHiggs 
code, we consider that all the proposed points are compatible with the LEP lower limit of 
113.5 GeV §. 

We also compute the rate for b s'-f decay and compare it with the experimental 



range 0. We implement the new NLO b ^ S'y calculations of |]T9| when M > 500 GeV, 
where M = mm{mq,mg). Otherwise, we use only the LO calculations and assign a larger 
theoretical error. For the experimental value, we combine the CLEO measurement with the 
recent BELLE resuh g], ^ 57) = (3.21 ±0.44 ±0.26) x 10"^. In our implementation, we 
allow CMSSM parameter choices that, after including the theoretical errors ath due to the 
scale and model dependences, may fall within the 95% confidence level range 2.33 x 10^^ < 
B{b 57) < 4.15 X 10"^ 

The final experimental contraint we consider is the — 2 value reported by the BNL E821 
experiment . This experiment has found an apparent discrepancy with the Standard Model 
prediction at the level of 2.6 a: 6a^ = (43±16) x 10^^*^. A large number of theoretical papers 



have discussed the interpretation of the BNL measurement within supersymmetry [^, ^1 



and they generally agree that /i > is favoured by the BNL measurement. The calculations 



we use in this paper are taken from [^, which are based on |2^, including also the leading 



two-loop electroweak correction factor p3 



We assume that R parity is conserved, and that the stable LSP is the lightest neutralino 
X p4|. We then constrain the CMSSM parameter space by requiring the calculated super- 
symmetric relic density to fall within the range 0.1 < ^^^"^ ^ ^■'^ preferred by astrophysics 
and cosmology. The upper limit on ^^^"^ conservative, being based only on the lower limit 
on the age of the Universe of 12 Gyr. Smaller values of Q^h"^ are certainly possible, since 
some of the cold dark matter might not consist of LSPs. However, allowing smaller values 
of Q^h'^ would open up only a very small extra region of the {mQ,mi/2) plane. 

Good overall consistency was found in between these relic density calculations, the 
LEP and other sparticle mass limits, the LEP Higgs limit, measurements of 6 — > 57 and the 
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recent BNL measurement of (7^ — 2, if > and tan/5 ^ 5. For tan/? > 50, there are not 
substantial regions with consistent electroweak vacua. 

The values of the CMSSM parameters for the benchmark points are shown in Table |I[ 
From these, soft masses are determined with ISASUGRA 7.51, and relic densities are cal- 
culated with a recent analysis |TT[ using SSARD that extends previous results [0] to larger 
tan (3 > 20. The chosen values of tan (3 range from 5 to about 50. In deference to (7^ — 2, most 
of the points proposed have fi > 0, but only about a half of the chosen points yield values 
of (7^ — 2 within 2a of the present central experimental value, and two of the points have 
/i < 0. The amount of CMSSM parameter fine-tuning required for electroweak symmetry 
breaking, along with the sensitivity of the relic density to the precise values of the input 
CMSSM parameters, are given in M together with the corresponding sparticle spectra. 



3 Direct Detection via Elastic Scattering 



The prospects for direct detection of neutralinos can be reduced to the computation of the 
neutralino-proton elastic scattering cross section. We first review the ingredients of this 
calculation that are implemented in SSARD. The MSSM Lagrangian leads to the following 
low-energy effective four-fermion Lagrangian suitable for describing elastic x-^ucleon scat- 



tering 25 



+ a^Xl^XQil^Qt + abiXXQil^Qi + a&iXl^XMi 



This Lagrangian is to be summed over the quark generations, and the subscript i labels up- 
type quarks (i = 1) and down-type quarks {i = 2). The terms with coefficients aii,a4i,a5i 
and a^i are velocity-dependent contributions and may be neglected for the purpose of di- 
rect detection calculations. The coefficients relevant for our discussion are, then, the spin- 
independent or scalar coefficients 
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and the spin-dependent coefficients 

1 
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4(mii - ml) 
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^Xal ~ 1^X4 1 



2 



4m| cos^ 

Here mu and are the squark mass eigenvalues, 



(3) 



2mwBi 



and the coefficients Z^. define the composition of the fightest neutrahno through 

X = Z^,B + Z^,W + Z^,H, + Z^,H2 . (5) 

The parameters e^, T^i, i/i denote electric charge, isospin and hypercharge (normalized so that 
ej = + ^), respectively, and 

Bi — (sin /3, cos /3) Ai — (cos /3, — sin /3) 

Cj = (sin a, cos a) — (cos a, — sin a) (6) 

for (up, down) type quarks. We denote by < the two scalar Higgs masses, and a 
is the Higgs mixing angle. Finally, rf^j^ are elements of the matrix that diagonalizes squark 
mass matrices through diag(m^j, m2j) = rj'' Mf {rf)~^ . 

The spin-independent (scalar) part of the cross section can be written as 

a^^^[ZU-r{A-Z)Uf , (7) 

TT 

where is the reduced neutralino mass and A, Z are the atomic number and nuclear electric 

charge, 

^= E fn^ + lMY.'^. (8) 

^ c "^o 27 frl 777.0 

P q=u,a,s q c,o,t y 

where 777p is the proton mass, and /„ is defined similarly. The parameters f^^ are defined by 

iT^pfrq = {P\mqqq\p) = rUgBg , (9) 
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while f^Q = 1 — J2q=u,d,s frq I13 . Following the analysis in [|^, we use the following 
values of f^^: 



cip) 



f: 



ip) 

Tu 



0.020 ±0.004 



Td 



0.026 ±0.005 



= 0.118 ±0.062 



where essentially all the error in f!^} arises from the uncertainty in strangeness composition 

25. 



(10) 



y 



Bd + Bu 

The corresponding values for the neutron are 



0.2 ±0.1 . 



(n) 
Tu 



An) 
JTs 



0.014 ±0.003 
0.118 ±0.062 



(n) 
Td 



0.036 ±0.008 



(12) 



These values are based in part on the experimental value of the 7r-nucleon a term 

a = ^(m, ± nid) x {Ba ± fi„) = 45 ± 8 MeV . 



(13) 



The larger value of cr = 65 MeV |^ considered by |^ leads to scattering cross section which 
are larger by a factor of about 3. It is clear already that the difference between the scalar 
parts of the cross sections for scattering off protons and neutrons must be rather small. 
The spin-dependent part of the elastic x-nucleus cross section can be written as 

32 



a2 = —G'pmiA'JiJ + 1) 



TT 



(14) 



where Gp is the Fermi constant, is again the reduced neutralino mass, J is the spin of 
the nucleus, and 

^ = -jMSp) + an{Sn)) , (15) 
where (Sp^n) are expectation values of the spin content of the nucleus and 



ip) 



r V2G, 



in) 



(16) 



The factors A^^'"''' parameterize the quark spin content of the nucleon. A recent global 
analysis of QCD sum rules for the gi structure functions ||32|, including 0{a^) corrections, 
corresponds formally to the values 



Ajf) = 0.78 ± 0.02 , A 



ip) 



-0.48 ±0.02 



A(f) 



-0.15 ±0.02 . 



(17) 
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In the case of the neutron, we have A^;') = A^"^ = A^f), and A^") = A^^). 

The calculation of the neutralino-nucleon elastic scattering cross-sections in Neutdriver 
is based on [^. The calculation of the spin-dependent contribution is identical to the one 
presented above. However, the spin-independent computation contains several additional 
pieces. First, the heavy flavor squark contribution is treated in exact one- loop approximation 
and (PD is replaced by 



as m 



P q=u,d,s H 



2 

Bd ^-^1^5 



where S^) and are loop integrals defined in Eqs. (18a) and (18c) of respectively. 
Second, Neutdriver also includes a Higgs contribution through squark loops (see the last 
term in Eq. (43) of |S^). Finally, it includes several contributions from twist-2 operators, 
which are listed in Eq. (46) of [Q. 

Numerical values from Neutdriver for the spin-independent and the spin-dependent 
components of the elastic cross sections for the scattering of neutralinos on protons and 
neutrons for each of the benchmark points are presented in Table |l|. (For other recent work 
in the CMSSM, see, e.g., 0.) In Fig. ||, we compare the results for the spin-independent 
(Tp and spin-dependent a'p cross-sections for neutralino-proton and neutralino-neutron scat- 
tering using SSARD |1^ and Neutdriver [14||. (For the latter, we have changed the default 
values of the quantities f^^, f^J, Af\ and Aj") to match those in (|10D, ([l|) and (^.) 
The differences are insignificant relative to the effects of different choices of CMSSM model 
parameters. Recall also that the mass spectra outputs of SSARD and ISASUGRA differ, as may 
be seen by comparing Tables 1 and 3 of |jl]]. Fig. |^ shows the projected sensitivities (a,b) 
for CDMS II and CRESST |B[ (solid) and GENIUS (dashed), and (c) a 100 kg 
NAIAD detector, as well as (d) the existing DAMA limit . Obtaining a competitive 
limit for the spin-dependent scattering on a neutron in the latter case might be possible with 
a large ^^Ge or Xenon detector. 

As was found in , there are strong cancellations in the spin-independent cross sections 
when fi < 0. These cancellations are due to sign differences between the up- and down-type 
quark contributions to the Higgs exchange terms in in (|^). Nominally, these cancellations 
occur only for a specific range in the neutralino mass. For tan (3 = 10, the cancellations occur 
for ~ 150 — 350 GeV, and are particularly effective when ~ 200 — 250 GeV. As one can 
see in Table |l| and Fig. |^, point D falls exactly into this range, thus explaining why its scalar 
cross section is anomalously small. Similarly, for tan (3 = 35, there are strong cancellations at 

~ 400 — 600 GeV |^ . Unfortunately, point K happens to fall in this range as well. Thus 



9 




Figure 2: Elastic cross sections for (a,h) spin-independent scattering and (c,d) spin- 
dependent scattering on (a,c) protons and (h,d) neutrons. The predictions 0/ SSARD (blue 
crosses) and Neutdriver (red circles) for neutralino-nucleon scattering are compared. Pro- 
jected sensitivities (a,h) for CDMS II and CRESST §^ (solid) and GENIUS /ff 
(dashed) and (c) for a 100 kg NAIAD array jj^/ , as well as (d) the existing DAM A limit 
are also shown. 



the two benchmark points with yU < are predicted to have very small spin-independent 
cross sections, but this would not generally be true for other CMSSM models with /i < 0. 
As one might expect, the differences between the SSARD and Neutdriver codes are largest 
for these points that exhibit delicate cancellations. 

Comparing the benchmark model predictions with the projected sensitivities, we see that 
spin-independent scattering seems to offer the best prospects for direct detection. Among the 
proposed benchmark points, models I, B, E, L, G, F, and C seem to offer the best detection 
prospects. In particular, the first four of these models would apparently be detectable with 
the proposed GENIUS detector. 
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4 Neutrinos from Annihilations in the Sun and Earth 



Dark matter particles collect in the gravitational wells at the centers of astrophysical bodies, 
leading to large densities and enhanced pair annihilation rates. While most annihilation 
products are immediately trapped or absorbed, neutrinos may propagate for long distances 
and be detected near the Earth's surface through their charged-current conversion to muons. 



High-energy muons produced by neutrinos from the centers of the Sun ||4l|] and Earth |^ 
are therefore prominent signals for indirect dark matter detection. 

The muon detection rate is dependent on both the neutralino annihilation rate and the 
resulting neutrino energy spectrum. The neutralino annihilation rate is proportional to the 
present dark matter density at the core of the Sun or Earth. Determinations of these densities 



are involved, but well understood. Various aspects of these calculations are reviewed in [14 



and estimates of neutralino annihilation rates in the CMSSM for both the Sun and the Earth 



are given in (For other recent work in the CMSSM, see, e.g.^ ||3, |4].) For the Sun, the 
annihilation rate has typically reached equilibrium and decreases for increasing neutralino 
mass. 

The neutrino energy spectrum depends on the neutralino composition. Neutralinos an- 
nihilate primarily to fermion pairs and gauge boson pairs. Annihilation to fermion pairs is 
helicity-suppressed, and so is significant only for heavy fermions, such as h quarks and r 
leptons, and t quarks if kinematically allowed. Neutrinos from these decays are typically 
rather soft. Annihilation to gauge bosons is possible only for neutralinos that are heavier 
than W bosons and have a significant Higgsino component. When possible, however, these 
annihilation channels typically dominate, producing hard neutrinos from two-body gauge 
boson decay. In this case, the muon flux is greatly enhanced, as both the cross section for 
conversion to muons and the muon range are proportional to the neutrino energy. 

Muon fluxes for each of the benchmark points are given in Fig. ^, using Neutdriver 
with a fixed constant local density po = 0.3 GeV/cm^ and neutralino velocity dispersion 
V = 270 km/s. For the points considered, rates from the Sun are far more promising than 
rates from the Earth. For the Sun, muon fluxes are for the most part anti-correlated with 
neutralino mass for the reason noted above. There are two strong exceptions, however: the 
focus point models E and F have anomalously large fluxes. In these cases, the dark matter's 
Higgsino content, though still small, is significant (see Table leading to annihilations to 
gauge boson pairs, hard neutrinos, and enhanced detection rates, as discussed above. 

The potentials of current and planned neutrino telescopes have been reviewed in ||15|| . 
The exact reach depends on the sahent features of a particular detector, e.g., its physical 
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Figure 3: Muon fluxes from neutrinos originating from relic annihilations inside (a) 
the Sun and (h) the Earth. Approximate sensitivities of near future neutrino telescopes 
= 10^ km"^ yr-i for AMANDA II NESTOR g^, and ANTARES g^, and 
$^ = 1 km^^ yr^"*^ for IceCube Ji4^J ) are also indicated. 



dimensions and muon energy threshold, and the expected characteristics of the signal, e.g., its 
angular dispersion, energy spectrum and source (Sun or Earth). Two sensitivities, which are 
roughly indicative of the potential of upcoming neutrino telescope experiments, are given in 
Fig. 1^. For focus-point model E, where the neutralino is both light and significantly different 
from pure Bino-like, detection in the near future at AMANDA II NESTOR |^, and 



ANTARES [0 is possible. Point F may be within reach of IceCube as the neutralino's 
significant Higgsino component compensates for its large mass. For point B, and possibly 
also points I, G, C, and L, the neutralino is nearly pure Bino, but is sufficiently light that 
detection at IceCube may also be possible. 

Muon energy thresholds specific to individual detectors have not been included. For 
AMANDA II and, especially, IceCube, these thresholds may be large, significantly suppress- 



ing the muon signal in models with less than about 4 to 6 E^^ ^ 



Note also that, 



for certain neutralino masses and properties, a population of dark matter particles in solar 
system orbits may boost the rates presented here by up to two orders of magnitude [^ . 
While this effect deserves further study, here we have conservatively neglected this possible 
enhancement. 
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5 Photons from Annihilations in the Galactic Center 



As with the centers of the Sun and Earth, the center of the galaxy may attract a significant 
overabundance of rehc dark matter particles |^|. Relic pair annihilation at the galactic cen- 
ter will then produce an excess of photons, which may be observed in gamma ray detectors. 
While monoenergetic signals from xx ~^ 77 and XX ~^ 7^ would be spectacular they 
are loop-suppressed and unobservable for these benchmark points. We therefore consider 
continuum photon signals here. 



%(E,.) = 5,6 X 10-» cm- x HdE^-^ ( '4) ( J(A(l) Afi . (19) 



The integrated photon flux above some photon energy threshold Eth is ||15| 

dNi (oiA dm GeV' ^ 

where the sum is over all annihilation channels dN^/dE is the differential gamma ray 
multiplicity for process i, AQ is the solid angle of the field of view of a given telescope, and 
J is a measure of the cuspiness of the galactic halo density profile. There is a great deal of 
uncertainty in J, with possible values in the range 3 to 10^ ||52|| . 

The integrated photon flux $(£'th) is given in Fig. ^ for each of the benchmark points. 
We choose AQ = 10"^ and a moderate value of J = 500. Estimates for point source flux 
sensitivities of several gamma ray detectors, both current and planned, are also shown. The 
space-based detectors EGRET, AMS/7 and GLAST can detect soft photons, but are limited 
in flux sensitivity by their small effective areas. Ground-based telescopes, such as MAGIC, 
HESS, CANGAROO and VERITAS, are much larger and so sensitive to lower fluxes, but are 
limited by higher energy thresholds. These sensitivities are not strictly valid for observations 
of the galactic center. Nevertheless, they provide rough guidelines for what sensitivities may 
be expected in coming years. For a discussion of these estimates, their derivation, and 
references to the original literature, see . 

Integrated fluxes for the benchmark points are given in Fig. ^ for two representative 
energy thresholds: 1 GeV, accessible to space-based detectors, and 50 GeV, characteristic of 
ground-based telescopes. Estimated sensitivities for two of the more promising experiments, 
GLAST and MAGIC are also shown. From (|1^), we expect the photon flux to be 
inversely correlated with neutralino mass. Roughly speaking, this general trend is seen in 
Fig. For Fig. Pd, it is offset by the requirement of a hard photon, which suppresses the 
signal from light neutralinos. In both cases, however, this general trend may be disrupted 
by a variety of additional effects. In particular, the photon spectrum is relatively hard for 
annihilation to gauge bosons; $(i?th) is, then, enhanced for the focus point models E and 
F, which have neutralinos with significant Higgsino components. The cross section cxj for 
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Figure 4: The integrated photon flux $(i?th) a function of photon energy threshold Eth 
for photons produced by relic annihilations in the galactic center. A moderate halo parameter 
J = 500 is assumed. Point source flux sensitivities for various gamma ray detectors are also 
shown. 



annihilation to bb through s-channel pseudoscalar Higgs is also enhanced for large tan/3 [Q, 
boosting photon signals at points I, J, K, L, and M. 

GLAST appears to be particularly promising, with points I and L giving observable 
signals. Recall, however, that all predicted fluxes scale linearly with J. For isothermal halo 
density profiles, the fluxes may be reduced by two orders of magnitude. On the other hand, 
for particularly cuspy halo models, such as those in all fluxes may be enhanced by two 
orders of magnitude, leading to detectable signals in GLAST for almost all points, and at 
MAGIC for the majority of benchmark points. 



6 Positrons from Annihilations in the Galactic Halo 

Relic neutralino annihilations in the galactic halo may also be detected through positron 
excesses in space-based and balloon experiments [0, |5^. The positron flux may be written 

as 



dVtdE 



(20) 



where p is the local neutralino mass density, the sum is over all annihilation channels i, and 
-Bg+ is the branching fraction to positrons in channel i. The initial positron energy distri- 
bution is given by the source function f{Eo), and the Green function G{Eo,E) propagates 
positrons in the galaxy. We use the Green function corresponding to a modified isothermal 



halo with size 4 kpc given in [^. The differential positron fluxes for the benchmark points 
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Figure 5: Comparisons between predicted integrated fluxes and prospective experimental sen- 
sitivities for photons with (a) a 1 GeV threshold, and (h) a 50 GeV threshold, following / (T^/ . 
Estimated sensitivities for (a) GLAST and (b) MAGIG /^]/ are also shown. A moderate 
halo parameter J = 500 is assumed. 



are given in Fig. |^. Note that the background spectrum drops rapidly with energy; hard 
positrons from neutrahno annihilation are most easily observed. 

To estimate the observability of a positron excess, we follow the procedure advocated 
in |]15|. For each benchmark spectrum, we find the positron energy ii^opt at which the positron 
signal to background ratio S/B is maximized. For detection, we then require that S/B at 
i^opt be above some value. The sensitivities of a variety of experiments have been estimated 
r5[| . Among these experiments, the most promising is AMS [|58|, the anti-matter detector 



m 



to be placed on the International Space Station. AMS will detect unprecedented numbers of 
positrons in a wide energy range. We estimate that a 1% excess in an fairly narrow energy 
bin, as is characteristic of the neutralino signal, will be statistically significant. 

Estimates of i^opt and the maximal S/B for each benchmark point are given in Fig. |^ 
To an excellent approximation, energetic positrons are produced only when neutralinos an- 
nihilate to gauge bosons that decay directly to positrons. Because this decay is two-body, 
Eopt ~ ""^x/^ foi' all benchmark points. As expected from (^0|) , S/B is typically inversely cor- 
related with neutralino mass. As discussed in Sec. ^for the case of neutrinos, however, there 
are two strong exceptions: the focus point models E and F. Again, these points have mixed 
gaugino-Higgsino dark matter. Rates for annihilation to gauge bosons and, consequently, 
the positron signals are therefore greatly enhanced. 

Even for points E and F, however, discovery of the positron excess is challenging for 
the smooth isothermal halo considered here. The positron search is most effective for light 
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Figure 6: Differential positron fluxes produced by relic annihilations in the galactic halo. 



Background fluxes are also shown for two models from iB^I^ 



neutralinos that are more Higgsino-like than those represented in this set of benchmark 
points. However, as with the photon signal, positron rates are sensitive to the halo model 
assumed; for clumpy halos |]59[ , the rate may be enhanced by orders of magnitude |]5^ . 



7 Conclusions 

In this paper, we have provided indicative estimates of the rates that could be expected for the 
benchmark supersymmetric scenarios proposed in |[^. We emphasize that, in addition to the 
supersymmetric model dependences of these calculations, there are important astrophysical 
uncertainties. These include the overall halo density, the possibility that it may be enhanced 
in the solar system, its cuspiness near the galactic centre, and its dumpiness elsewhere. 
For these reasons, our conclusions about the relative ease with which different models may 
be detected using the same signature may be more reliable than the absolute strengths we 
predict, or comparisons between the observabilities of different signatures. Nevertheless, our 
estimates do indicate that there may be good prospects for astrophysical detection of quite 
a large number of the benchmark scenarios. 

In particular, the direct detection of relic particles by spin-independent elastic scatter- 



ing in models I, B, E and L may be possible using the projected GENIUS detector, 
with models G, F and C not far from the likely threshold of detect ability. The prospects of 
detecting spin-dependent elastic scattering do not, however, look so promising in the bench- 
mark scenarios studied. The indirect detection of muons generated by high-energy neutrinos 
due to annihilations inside the Sun should be most easily detectable in models E, F and 
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Figure 7: (a) The optimal energies for which the positron signal to background ratios S/ B are 
maximized, and (h) S/B at these energies for each of the benchmark points, following /(J5|/. In 
(a), the dashed line is for E^pt = rn^/2, and in (b), the estimated sensitivity of the AMS 
experiment is shown. 



B, followed by models I, G, L and C, which offer prospects with the proposed IceCube 
detector. However, unless there is a substantial solar-system enhancement, the prospects for 
detecting annihilations inside the Earth are not so encouraging. Models L and I offer the 
best prospects for the detection of photons from annihilations in the galactic centre, followed 



by models K, B, E and G. Here the best prospects may be those for the GLAST satellite, 
with its relatively low threshold. However, there may also be prospects for ground-based 



experiments such as MAGIC |54|, if the halo is cuspier at the galactic centre than we have 
assumed. Models E, F and B offer the positron signals with the largest signal-to-background 



ratios, though apparently requiring a sensitivity greater than that expected for AMS 
unless the halo is rather clumpy. 

In specifying the benchmark models, the constraint coming from the anomalous magnetic 
moment of the muon was not imposed rigorously. However, it was noted that the more 

— 2-friendly models I, L, B, G, C and J offered good prospects for detecting several 
super symmetric particles at the LHC and/or a linear e^e~ collider with 1 TeV in the centre 
of mass. Most of these models also exhibit good prospects for dark matter detection, with 
the exception of model J. Among the less g^ — 2-friendly models, we note that E, F and 
K offer some astrophysical prospects. This is particularly interesting in the case of focus 
point model F, which does not offer generous prospects at colliders, and model K, which 
is not easy to explore with a linear e^e~ collider. On the other hand, models M and H, 
which are particularly difficult to explore with colliders, also do not offer bright prospects 
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for astrophysical detection. 

Our analysis indicates the effort required to cover the possible supersymmetric param- 
eter space via a number of different astrophysical signatures, at least within the CMSSM 
assumptions used here. It would be interesting to extend such a benchmark analysis to other 
types of supersymmetric models, but that lies beyond the scope of this paper. Ultimately, 
one would hope to be able to confront accelerator and astrophysical measurements of super- 
symmetry, and make non-trivial cross-checks of our CMSSM assumptions, but that is for the 
future. For the moment, the race to discover supersymmetry is still open, and our analysis 
indicates that there may be good prospects for detecting supersymmetric dark matter before 
the LHC comes into operation. 
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